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TECHNICAL NOTE 2310 

GENERALIZATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS TO 
THREE-DIMENSIONAL FLOWS INCLUDING THOSE OF ROTATING SYSTEM 

.By Artur Mager 


SUMMARY • 

The Navier-Stok.es equations of motion and the equation of continu- 
ity are transformed so as to apply to an orthogonal curvilinear coordi- 
nate system rotating with a uniform angular velocity about an arbitrary 
axis in space. A usual simplification of these equations as consistent 
with the accepted Boundary- layer theory and an integration of these 
equations through the Boundary layer result in Boundary-layer momentum- 
integral equations for three-dimensional flows that are applicable to 
either rotating or nonrotating fluid Boundaries. 

These equations are simplified and an approximate solution in 
closed integral form is obtained for a generalized Boundary-layer 
momentum-loss thickness and flow deflection at the wall in the. turbu- 
lent case. 

A numerical evaluation of this solution carried out for data 
obtained in a curving nonrotating duct shows a fair quantitative agree- 
ment with the measured values. 

The form in which the equations are presented is readily adaptable 
to cases of steady, three-dimensional, incompressible Boundary-layer 
flow like that over curved ducts or yawed wings and it also may Be used 
to. describe the Boundary-layer flow over various rotating surfaces, thus 
applying to turbomachinery, propellers, and helicopter Blades. 


INTRODUCTION 

The development of the Boundary layer on the various parts of 
turbomachinery (compressors and turbines), helicopter Blades, propel- 
lers, and in curved ducts is influenced By centrifugal and Coriolis 
forces in addition to the pressure and viscous forces. As a result of 
these forces, the flow in the Boundary layer not only has the char- 
acteristic velocity deficiency But also has. Because of this velocity 
deficiency, direction different from that of the flow outside the 
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boundary layer. Thus the behavior of the boundary layer in three- 
dimensional flow may be quite unlike the behavior in two-dimensional 
flow. The effect of these additional forces on the boundary layer has 
been realized for some time and the observed discrepancies in the 
boundary-layer behavior have usually been explained only in a qualita- 
tive mann er as, for instance, in references 1 to 4. 

The literature concerning the theoretical aspect of the three- 
dimensional boundary-layer flow is meager. For the laminar case most 
of the published work has been carried out in connection with the yawed 
wing (references 5 to 7). For the turbulent case, although a number of 
researchers have established the general form of the differential equa- . 
tions applicable, no actual solutions of these equations have been 
obtained. Tetervin, for instance, presents boundary-layer momentum- 
integral equations in three dimensions for a fluid of variable density 
and viscosity (reference 8). Gruschwitz establishes -the momentum- 
integral equations for boundary- layer flow along an arbitrarily curved 
streamline in reference 9. Burgers gives the differential equations, 
on the development of boundary layers in the case of axially symmetric 
flows having a rotational component (reference 10) . Prandtl, in addi- 
tion to presenting a form of three-dimensional momentum- integral equa- 
tions, suggests the general procedure that could be followed to obtain 
a solution (reference ll) . Experimental data are similarily lacking. 

In spite of considerable literature search, only the data of Gruschwitz 
(reference 9) for a curved duct and the data of Kuethe, McKee, and 
Curry (reference 12) for a yawed wing were found. 

As a result of research on this problem at the NACA Lewis labora- 
tory, the boundary- layer momentum- integral equations are derived and 
presented herein for a set of orthogonal curvilinear coordinates, which 
may or may not be rotating about an arbitrary axis in space and can be 
laid out along a streamline of the potential flow. The so generalized 
equations are then transformed by use of an assumed velocity distribu- 
tion and friction law for turbulent boundary layer so that an approximate 
solution can be obtained for the boundary-layer momentum thickness and 
the direction of boundary-layer flow. Finally, a numerical solution is 
carried out for the Gruschwitz data in order to make a comparison 
between the estimated and actual measured values. 

The equations as given in their generalized form are readily 
adaptable to cases of steady, three-dimensional, incompressible boundary- 
layer flow, involving centrifugal and Coriolis forces. The approximate 
solution, however, has been carried out only for the turbulent boundary 
layer, because in most of the aerodynamic configurations, where these 
equations apply, transition from laminar to turbulent flow occurs 
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comparatively early in the flow process. A laminar form of the approx- 
imate solution can he obtained by simple substitution of a suitable 
velocity profile and friction law. 

It should be noted that whereas the differential equations describe 
the flow phenomena with only the accepted simplifications, the approxi- 
mate solution depends to some extent on the assumed boundary- layer 
velocity profiles and the relation for friction. Both of these assump- 
tions were made on the basis of the data of Gruschwitz (reference 9) 
only, because the data of reference 12 were not adaptable to extensive 
computations for the purpose of this analysis. The measurements of 
Gruschwitz, on the other hand, have certain shortcomings as they were 
obtained in a nonrotating channel formed by two circular-arc shaped 
walls . Thus the generality of the velocity profiles measured by 
Gruschwitz is in question. A revision of the approximate solution can 
therefore be expected when more data become available. In addition, 
any speculation on the occurrence of boundary- layer separation (which 
by definition is a special form of a velocity profile) would be abso- 
lutely meaningless j no further mention will therefore be made of this 
phase of the problem. 


SYMBOLS 

The following symbols are used in this report (the dimensions are 
given in right-hand column) : 


A 

constant occurring in second approximation for 0 

O' 1 ) 

a 

resultant acceleration vector in fixed (inertial) 
system 

(Zf 2 ) 

*o 

resultant acceleration vector in Cartesian coor- 
dinate system 

(It" 2 ) 

B 

constant occurring in second approximation 
for 0 

(o) 

b 

position vector of particle 

0) 

c 

curvature of x-axis (fig. l), ^ 

O' 1 ) 

d 

constant > & 

0) 
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E-l(x) 

E 2 (x) 


e 3 U) 


e xx' e xy 


F 

f X> f Y> f Z 

G 


H,J,K,L 

l 

m 


m 

P 

P XX' P YY' * ’ * 


0.01255 

(K-J) 


dx 

0 T 

xi I, 


£ / c€ dx ' 

^ /xi 


18 T 


*i 


°v 

C ^ dx, 


rate-of -strain components 

resultant-force vector acting on particle 
components of tody forces per unit mass 


( 0 ) 

• ( 0 ) 

( 0 ) 

(f 1 ) 

(mlt" 2 ) 

(It" 2 ) 


function describing boundary-layer velocity profile, 
also taken as ( 


( 0 ). 


function describing boundary-layer velocity profile, 

also taken as ^1 - ^ (0) 

quantities describing relations among various char- 
acteristic loss thicknesses in boundary layer (0) 

transformation coefficients (o) 


length 

parameter determining nature of boundary- layer 
equations 

mass ' 

static pressure 


( 0 ) 


(ml ^t 2 ) 


components of stress per unit area in Cartesian _ 2 

coordinate system (ml t ) 
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% 

R . 

R(e x ) 


s 


s 

t 


■U,v,w 


U 0 ,Vo,W o 

u,v,w 


X,Y,Z 


x,y,z 

z x (x) 

a 


p 

6 


& a c -JT 

X > & Z 


£p'->£+>£- 


resultant velocity vector 

perpendicular distance of particle from axis of 
rotation 

e x u 

Reynolds number based on ®x> IT 


radius of circle 


total path length 



arc length 
time 

values of u, v, and w outside boundary 
layer 

/ 

velocities in Cartesian coordinate system 


(it" 1 ) 

(i) 

( 0 ) 

(i) 

(i) 

(i) 


(Zt-i) 

(It -1 ) 


time averaged velocities in curvilinear coordinate 
system 

Cartesian coordinate system 
orthogonal curvilinear coordinate system 
function used in transformation 


(if 1 ) 

(i) 

(i) 

(i) 


boundary-layer deflection angle measured from 
direction of resultant skin-friction stress to 
direction of flow outside boundary layer (o) 

angle between X-axis and tangent to x-axis (0) 

boundary- layer thickness - (z) 

displacement thicknesses in three-dimensional 

boundary layer (z) 

measure of boundary-layer deflection, tan a (o) 

slope of characteristic line (o) 
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0 ' 

X 

V 

p 

°x-’ T xy ’ ' * 


generalized boundary-layer momentum-loss thickness, 

e x R 1 / 4 (0x) ' (*) 


N 

momentum- loss thicknesses in three-dimensional 
boundary layer 

/ ,^ x 

variable of function i|r. 


r x 

S-J / 

J x i 


c dx , 


(0 

(0 


kinematic viscosity 
components of vorticity vector 
density 

apparent stresses existing in turbulent flow 


(zV 1 ) 

(f 1 ) 

(ml- 3 ) 

(mZ“^t - ^) 


T . shear stress at wall (mZ _ ^t _ ^) 

o'. 


t 


U3 


CO ,03 ,03 

x' y z 


Subscripts : 


arbitrary function satisfying 
equation (39) and boundary 
conditions 



angular velocity 

components of vector 
system . 

components of vector 
system 


(f* 1 ) 

co ‘ in Cartesian coordinate 

(t- 1 ) 

co in curvilinear coordinate 

(t- 1 ) 


i initial value 

x x-direction 

z z-direction 

IjII order of approximations 

For Gruschwitz data-point designations and streamline designa- 
tions, see figure 2. 
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DERIVATION OF BOUNDARY-LAYER MOMENTUM- INTEGRAL EQUATIONS 


The equations for steady flow of a fluid having constant density 
are derived in a Cartesian coordinate system X,Y,Z, rotating with 
uniform angular velocity about an arbitrary axis in space. These equa- 
tions are then transformed to an orthogonal curvilinear coordinate 
system x,y,z, such that the x-axis can be placed along any convenient 
path in the XY -plane, which is considered as a plane of a wall. These 
equations are then simplified in a manner consistent with the boundary- 
layer theory. If. the path is chosen so as to match a streamline of the 
potential flow, only one velocity will exist outside the boundary 
layer, that along the streamline. Furthermore, the changes in boundary- 
layer quantities in a direction other than that along the streamline 
are expected to be relatively small in comparison to the changes along 
the streamline. Additional simplifications may thus be possible. 
Finally, integration through the boundary layer gives the generalized 
form of momentum- integral equations for three-dimensional flows that 
may or may not involve rotation of the system. 


Equations for steady flow of fluid with constant density in rotat- 
ing Cartesian coordinate system. - The Navier-Stokes equations of flow 
for a fixed Cartesian coordinate system X,Y,Z (reference 13, p. 576) are 


DU 0 

Dt 

DV 0 

Dt 

DW 


dp 


XX 


dpyx dp 


ZX 


- p f x + -ar + ~w + ~5z 

dpxY dpyy dp ZY 

- p f Y + ~ar + + “Sz - 


5p xz , forc , dp zz 


P DV = pfz + ax 


5y 


az 


(la) 

(lb) 

(lc) 


and the equation of continuity is 


au Q av G aw 0 

w + w + az~ = 0 


( 2 ) 


It is now assumed that this Cartesian coordinate system is rotating 
with a uniform angular velocity co and that the observations of the 
motion of the fluid particles are still made_from a position rigidly 
attached to the same system. The velocity q 0 and acceleration a Q 
are as seen by the observer, that is, they are relative to the X,Y,Z 
system. Because of the rotation, however, the X,Y,Z system is not an ' 
inertial system (reference 14, p. 53) and_thus the second law of motion 
holds only with respect to acceleration a relative to some other 
system that is nonrotating. 
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ma = F 

In terms of a Q then (reference 14, p. 104), 

ma 0 + men X (u) X b) + 2mn X = F 

Here mo X (o> X b) represents the centrifugal force and 
2mm x^ is the Coriolis force. 


Thus for a Cartesian coordinate system rotating with a uniform 

DU q DVq • DW 0 

angular velocity <n, the expressions for and must 

be modified by proper components of the Coriolis and centrifugal accel- 
erations. For steady flow, the component accelerations as referred to 
a rotating Cartesian coordinate system are therefore 


DU 0 

Uo 

Su 0 


Vo 

a 

o . 


Wo 

Dt 

Sx 

+ 

Sy 

+ 

RV o 

U o 

SV o 


v o 

SV o 


w 0 

Dt _ 

Sx 

+ 

sr 

+ 

DW 0 

u o 

Sw 0 


v 0 

Sw G 


w o 

Dt 

Sx 

+ 

sr 

+ 


UUq 2 

+ 2 ( a YWo-<*>zV 0 )-a> R 
Sv Q 

3zT + 2 (“z U o-° i K W o)- a)2R 

sw D „ 

+ 2 ( a >x v o -a> f U o)" a) R 


Sr 

Sx 


Sr 

Sy 

Sr 

Sz 


(3a) 

(3b) 

(3c) 


The equation of continuity, which does not involve any accelerations, 
remains the same. 


Transformation to orthogonal curvilinear coordinate system. - 
Transformations similar to those of Gruschwitz (reference 9) are used 
as in dicated in figure 1 with the precaution that the system remain 
right-handed. 

f X 

X = I cos p dx + z sin p 

JO . 

x = y > (4) 


where 


'X 


Z = Zi + z cos p 


J 


Z-j_ = constant - 



sin p dx 
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and 

P » POO 

Use of these transformations permits an arbitrary curvature of the 
x-axis in only one plane, the XZ-plane. Thus the solution is somewhat 
restricted. In two-dimensional boundary-layer investigations, however, 
it is found that the boundary-layer equations are unaffected if the 
radius of curvature in the XY-plane is large as compared with the 
boundary- layer thickness (reference 15, p. 120). In three-dimensional 
boundary layer the same limitation will probably apply providing, of ' 
course, the values of ay ay . and a) z are properly adjusted to take 

care of this additional curvature. Setting 



(curvature of x-axis) 



Thus, 

(ds) 2 = (h 1 ) 2 (dx) 2 +(h 2 ) 2 (dy) 2 +(h 3 ) 2 (dz) 2 =(<3X) 2 + (dY) 2 +(dZ) 2 
But because 

dX = ^dx+^dy + ^dz 

' and so forth,. 

(ds) 2 = (l+cz) 2 (dx) 2 + (dy) 2 + (dz) 2 
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and 


h-j_ = (l+cz) 


h 2 = . 1 


h 5 = l 


( 5 ) 


The expressions for the linear accelerations can be written directly, 
as given in reference 13 (p. 158). (it should be noted that the h 

values herein are reciprocals of those in reference 13.) -The components 
of a gradient now are 


_i_ a 

h^_ 3x 


_i a 
h 2 3y 


i_ a 

hj S? 


whereas the components of ^ X <n remain 


Vcd„ 


-WCDy 


Wo^-Uo^ UOJy-VflDjj 


Thus the accelerations in the rotating x,y,z system are written as 


du u au v au w au 


V 


ah n 


ahc 


w 


ahn 


ah. 


ST - v *r) + sr - » arV + 

And the expressions for ^ and ^ follow from symmetry. The equa- 
tion for the divergence now has the form 

aiT 5o “ [s (n 2 b 3 u) + (hjijV) + £ (hiVOJ - o 

Whereas the components of the curl ^ are 

* = h^ (h 3 W) - ^ (h 2 V 3 
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In order to obtain the viscous terms the preceding expressions axe used 
in the expansion of 


v^grad (div q^ - curl (curl q^J 


If equations (5) are substituted into these general expres- 
sions and the differentiations are carried out, the equations for flow 
in orthogonal curvilinear coordinate system rotating with an angular 
velocity co are obtained. The body forces are neglected here. 

Tfb H + v 1f + w H + m ~ ikz cu2r H + 2 (v^z v) 


1 

1 dP 

- 4 - 

r i 

d 2 U z 

dU dc 

d 2 U 

1+cz 

p dx 

* 17 

[_(l+cz)? 

dx 2 (l+cz) 

3 dx dx 

dy 2 

d^ 

c 

dU 

Uc 2 

W dc 

2c 

dW 

dz 2 

1+cz 

dz 

(l+cz)^ 

(1+cz) 3 dx 

+ 

(l+cz 

)2 dx 


(6a) 


U 


dv 

dx 


„ dV tI dV 1 2 d dR _/ TT T7 x 1 dP 
+ v 57 + w 51 - I+SS “ R d^ + 2 K U -x W > = " p dy + 


dV 


dR 


1 dP 


d 2 V 


_ . z dV dc . d^V d^V c dV 

|_(l+cz) 2 dx 2 " (l+cz) 3 ^ dx dy 2 + dz 2 + 1+cz ^ 


] 


(6b) 


U dW . „ dW dW c „2 2^ dR o/ „ lA _ 1 dP . 

1+zc dx dy dz 1+zc dz j ' p dz 


d^W 


U dc 


1 dU dV dW Wc 

1+cz 3x <5y 3z; 1+cz 


z dW dc 


•cz) 3 Six . 


2 c dU~l 

(1+cz) 2 

(6c) 

= 0 

(6d) 


In the general orthogonal coordinates, the expressions for the rate-of- 
strain components axe 
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1 Su _v_ Sh i w_ dh i\ = ^3 a /_w_\ + l£ a t-L\ 

Gxx ~ \h^ Bx h^hg Sy ^ 3^1 az / h 2 \h 3 / ' h 3 3z l^hg/ 


Sh n 


e = 

yy 


'ZZ 


- zi k 


‘ 2 (si i + 


dV W Sh 2 U ^ h 2 

+ " — r + 


0V + w £ + - - 
3y h 2 h 3 dz h-jh 2 > 

U ^ h 3 + V ^ h 3\ 

h 3 h! 3T - + h 2 h 3 3y~/ 


X 1 S /U \ . n 3 b / W 


e zx h 3 3z ^h ly / hi 5x \h 3 

= ^ a (J_\ h. b /U_ 

e xy h-L 3x lh 2 y h 2 5y \h x 


The viscous terms in equations (6a ) } (6b ) } and (6c) may be expressed 
using the rate-of-strain components as 


V 




&J dc . d 1 2 U . a 2 U 


au 


Ucr 


( 1+c z ) 2 Bx 2 " (l+cz) 3 Sxdx ay 2 ^ z 2 l+cz 5z (l+cz)‘ 


W 


(l+cz) 


_ dc 
3 dx 




de T 


= V 


(l+cz) 

2 3x_ 

= V 

[l+cz dx 

, */- 
+ _ ^7 

a 2 v 

z 

av dc 

a 2 v . a 2 v 

ax 2 

(l+cz) 3 dx 

ay 2 

az 2 

^ e yy 

1 

ae X y 

de yz 

c 

_ay 

l+cz dx 

~5T + 

l+cz ' 


^ e zx 2ce xz 


bV 


(7a) 


(7b) 


Wc‘ 


1 b^Vt U dc z 5w dc + c) 2 W ^ b^W wu ^ 

” fp^j 2 5? ' (l+cz) 3 S ‘ (l+cz) 3 S S 3y 2 3z 2 ‘ (l+cz) 2 


c aw 

l+cz ciz 


2c _ au 

2 chc 


(l+cz) 


■1= v[li 

: [l+c 


- ae xz. ae yz 5e zz 
l+cz ax • ay az 


l+cz ( e xx _ e zz) 


(7c) 


Equations (6) are directly applicable to the laminar flow. For tur- 
bulent flow, because of the velocity fluctuations it is necessary to 
modify the stresses by addition of the so-called Reynold's stresses. 
Thus, maxing use of the parallel form in equations (7), the Navier- 
Stokes equations of motion for turbulent flow may be written in terms 
of the apparent stresses as 
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_U_ ^ + + wS + r-2— UW. - 

1+cz dx oy oz 1+cz 


A, M 


— - « || + 2(o)yW-(O z V) 


crfR 


1 i 3 p 


_ _ I 1 ^°x ^ T yx ^ T zx 

1+cz p 3x + p |_l+cz dx + dy + dz + 1+cz 

u I + v I + “ I - ub‘“ 2R I + 

1 dP j_ 1 r 1 ^ T xy 

“ - 3TT ~ 111 


2c T zx 
1+CZ _ 


den, 


dr„ 




(8a) 


(8b) 


U dW 
1+cz 3x 



+ W 


cJW c 
3z ” 1+cz 


U 2 -<jd 2 R 


3£ + 2(c0 x v ^Y u) 


1 dP 1 1 ^xz ^ T yz + ^z _ c 

p 3z + p _l+cz dx dy dz 1+cz 



(8c) 


Simplification for flow within boundary layer. - Equations (6) and. 
(8) are equivalent to the complete Navier -Stokes equations . Within the 
bou ndar y layer , however, certain terms, whose contribution is relatively 
unimportant can be neglected. If the y-axis is taken as normal to the 
wall, the boundary-layer flow then takes place over the xz-plane (or 
the XZ-plane) . All, terms are now made dimensionless by referring the 
lengths to some body length, the velocities to their free -stream 
values, an d so forth, as explained in reference 16 (p. 45) and all 
quantities of the order of magnitude of 5 or smaller are neglected. 
Furthermore, because the boundary-layer flow along a definite path 
z =0- is of interest, additional simplifications are possible. Setting 
z = 0 restricts the equations, because the general boundary conditions 
(not on the x-axis) cannot be satisfied. It will subsequently be seen, 
however, that these general boundary conditions are unnecessary in the 
solution of the final equations. These simplifications yield the 
Navier -Stokes equations for flow within the boundary layer in a 
rotating orthogonal curvilinear coordinate system evaluated at z = 0, 


u 


du 


+ v 


du 

5y 


+ w 


du 

31 


cuw - - — 
p 


dP 

3£ 


+ CD 


^ i - ^ + v (0) <9a) 


- uy 


2 R || +. 2 (cD 2 U - CDjjW) 


dP 

3y 


(9b) 
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u 


dw , 


dw 


dw 


- u 2 c = 


dp 

3I + 


0) 2 R 


3r 

SI + 


Zoyi 


+ v 



(9c) 


for the laminar case. For the turbulent boundary layer, a corresponding 

1 oT X y 

set of equations, is obtained with the substitution of -r-^ 

, dT /I o \ p °y 

, 1 z y 

and — !■— — 


for 


■&) 


P oy 


for 


•(&)■ 


Equation (9b) shows, as pointed out in reference 10, that because 
all the terms on the left-hand side of the equations are of the order of 
magnitude of one, within the boundary layer, P can vary at most by an 
amount of the order of 6. It is reasonable then to neglect this varia- 
tion and consider P solely, a function of the flow outside the boundary 
layer. Thus, if x is chosen to coincide with a streamline of the flow 
outside the boundary layer, V = W = 0, and by integration of equa- 
tion (9a) with the effect of viscosity neglected the following relation 
is obtained: 

P = constant - ^ pU 2 + ^ pao 2 R 2 (10) 


which is a form of the equation of Bernoulli. 


Furthermore, because outside the boundary layer the flow with 
respect to some nonrotating set of coordinates is irrotational, with 
reference to the rotating coordinates the components of the vorticity 
vector become 


| = -2cu x 
T) = -2u>y 
t “ -2®z 


(11) 


This assumption of irrotationality is not always true and in some 
applications, such as the later stages of an axial compressor, it cannot 
be used. As long as vorticity is distributed according to some definite 
pattern, however, a relation between the components of vorticity and the 
components of rotational velocity may be found and substituted for 
equations ( 11) . 


Substituting again in the expression for the components of vor- 
ticity gives 

= (1+cz) u • s"] 
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And for z = 0, W = 0,' 
for curvature becomes 


which is along the streamline, the expression 


c ss 



1 dtA 

U Sz / 


( 12 ) 


The equation of continuity remains 


du , dv . 

+ 37 + 



wc 


0 


(13) 


Generalized boundary-layer momentum-integral equations. - In order 
to obtain the boundary-layer momentum- integral equations, equations (9a) 
and (9c) are integrated with respect to y through the boundary layer 
to some constant height d such that 


d > 6 


a 

3x 


r d 

u 2 dy + 2a>y 
jO u 


wdy + |! 

0 dz J 


w dy - U I u dy - U / w dy + 


' 0 


ai 




a 




Pd 


uw dy - 4 


u 


, „ i au 

uw dy - 2 - ^ 


i ^ au T c 

uw dy + 2cDy W dy = Ud ^ 

10 

(14a) 


and 


a 


Pd 


uw dy + 


0 




Pd 


w 2 dy + c 
0 u 

^d 


w^dy - c 1 u^dy 


0 


au 


= U ^ dy + 2ufy 


, °t Z 

u dy ■ —£~ 


(14c) 


These equations apply equally well for the laminar or turbulent boundary 
layer, with the value of T Q representing the shear stress at the wall 
accordingly adjusted. By suitable use of equations (12) and (13), these 
equations may be transformed to 
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JL (U-u)u dy + gj o (U-u)dy * £ J q (U-u)w dy 


d Pd 

4w T 0 

(U-u)w dy rf (U-u)w dy = — 


u^J 0 - U J Q 


P 


(15a) 


s/o ” 47 + ®Jo '** + t sjo " 2iy . 

■ u dy + **fo u "" *'%fo “ 2dy ' *J. uZdy ' ^ 


The following definitions are now introduced: The momentum 

thickness in the x-direction of the flow in the x-direction, 


(15c) 


6 =— (U-u)udy 

u n 


The displacement thickness in the x-direction, 


5* ■ jj I (U-u) dy 


momentum thickness in the z-direction of the flow in the z-direction, 


6 Z = w 2 dy 

IT „ 

J 0 

The displacement thickness in the z-direction, 


&z = u W dy 
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The momentum thickness in the z-direction of the flow in the x-direction, 

r d 

(U-u)w dy (20) 

0 


e xz - ^ 


The momentum thickness in the x-direction of the flow in the z-direction, 

'd 


Q z x = -p I m d y 
IT 


( 21 ) 


All these thicknesses, as in two-dimensional boundary-layer 
theory, have a dimension of length. Furthermore, 


5* - 0 Y17 = i | w dy - i I (U-u)w dy = — j uw dy = 6 


xz u 


U 2 l 0 


U 2 . IQ 


7 ZX 


(22) 


With the use of definitions (16) to (21) and equation (22), equa- 
tions (15a) and (15c) reduce for z = 0, to 


S0 x 1 dU w 

+ U (2e x + 5 x ) + 


^0 CO ' *r 

_2E. _ 4 JL e _ 

* u xz „2 


pU 


(23a) 


and 


^®Z "®xz), i 

‘ ^ "T” TT 


dz 


ox 


o i < > * ! i ^ <»*+«*> - - T -¥ 

P U 

(23c) 


Reduction of equations to forms obtained by other investigators . - 
If only two-dimensional flow exists, that is if c = 0, w. = 0, and 
00 =, 0, then equation (23c) vanishes and equation (23a) becomes an 
ordinary Karman momentum-integral equation 

pU 

If (D = 0, that is, if the system is nonrotating, equations (23a) 
and (23c) become identical with the equations of Gruschwitz 
(reference 9) . 
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Setting c = 0 in equations (9a) and (9c) makes these equations 
identical with the equations of Burgers (reference 10), who carried 
out his derivation for a Cartesian coordinate system. 

Finally, if the system of coordinates is chosen so as to maintain 

the right-hand rule and c is set equal to - — , thus establishing the 

dz r 

x-axis as a circle, then ^ = 1 and because of axial symmetry all 

derivatives with respect to x vanish. The coordinates are now 
assumed to be in a fluid that is motionless at great distance from the 
surface of the rotating immersed disk. Thus, 


U = W = 0 


CD = 0 


Integration of equations (9a) and (9c) gives after some manipu- 
lation, 

2 * If (* 2 J uv dy ) = - ^ 2 * r2 

\ 

and 


I w2dy )- I u 2 ay- --2^r 


'0 


0 


which are identical with equations of von Karman for the rotating disk 
(reference 17). 


APPROXIMATE SOLUTION OF MOMENTUM- INTEGRAL EQUATIONS 


FOR TURBULENT BOUNDARY LAYER 

Transformation and reduction of dependent variables. - In order to 
obtain a solution of. the momentum- integral equations, additional rela- 
tions are needed describing the velocity profiles existing in the 
boundary layer and the friction at the wall.- 

With the use of a suggestion by Prandtl (reference 11), the 
expressions for u and w that will be used are 
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u = UG 
w = £UG 


(I) 

m) - 


eug 


(24) 


with boundary conditions 

on 

G 

and g 


for 

y = 

5, 

G = 1, 

g = o 

for 

y ® 

o. 

G = 0, 

g = 1 

t 

and with € defined as 







€ 

= tan a 



(25) 


where a is the angle between the direction of the resultant skin- 
friction stress and the direction of the flow outside the boundary 
layer. Because of this definition of € , g-1 at y=0 because 


c)w 

lim e 

ou 

y~>° ^ 


or 


T 


o, z 


€ T 


o,x 


(26) 


Mathematically,' such use of € implies a linear variation of w‘ 
with e arid makes possible the dissociation of the w velocity profile 
from its scale and direction. Because the flow must change direction 
in the boundary layer from that at the wall to that in the free stream, 
there is no reason to assume that such a dissociation is actually pos- 
sible. In other words, there is no reason to believe that g should be 

a function of (y/s) alone and not of £ as well. In accordance with 

reference 11, however, this approx im ation is certainly admissible for 
small values of 6 and gives results of qualitative accuracy for 
moderately large € . In addition, in order to check this assumption, 

the value of ' i Z, for several experimental velocity profiles and 

values of € ranging from 0.216 to 0.670, as obtained from refer- 
ence 9, is plotted against y/& in figure 3. The results of this plot 

indicate indeed that Gg is independent of e. . 
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In parallel to the two-dimensional boundary-layer theory, the 
following definitions are made: 

'd 

(1-G) dy 

0 

= H 


0 


(l-G)G dy 


(l-G)Gg dy 


0 


U-G)G dy 


s J 


Gg dy 


(l-G)G.dy 


= K 


G 2 g 2 dy 


0 


(l-G)G dy 

0 


= L 


\ 


(27) 


The relations among the various thicknesses may then be written 

* ^ 

6* =H0 X 

0 xz = £ J 0 x 


S* = e K 0 x 


e z = € Le x 


(28) 


J 


The quantities H, J, K, and L are functions of G and g. 
Because G and g -are representative of the velocity profiles in the 
boundary layer, the changes in these velocity profiles must be reflected 
in turn in the values of H, K, J, and L. In other words, the 
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external forces acting on the boundary layer and influencing the changes 
in the shape of the velocity profiles also cause a variation in ' H, K, 
J, and L. Unfortunately, the available data of reference 9, do not 
involve large changes in the shape of the velocity profiles and the 
quantities H, J, K, and L. This behavior of the velocity profiles 
is verified in figure 3. The data of reference 12 do indicate large 
changes in the shape of the velocity profile j however, the data are 
not presented with sufficient detail to permit an accurate evaluation 
of H, J, K, and L. Thus, until more extensive experimental data 
become available, the quantities H, J, K, and L are assumed to be 
constants that can be evaluated either by assuming a suitable form for 
G and g or by computing directly from Gruschwitz data. 


In accordance with reference / 9/ good assumptions for G and' g 

are: 



(29) 


An indication of the degree of fit afforded by these expressions 
can be obtained from figure 4(a), -where a calculated profile with 
n = 7 is compared with one of the profiles of Gruschwitz. Other pro- 
files of Gruschwitz data give similar results. It should be noted that 
this good agreement should not be interpreted as meaning that assump- 
tions (29) will always give a good representation of the velocity pro- 
files in the three-dimensional turbulent boundary layer. Figure 4(b) 
shows. a comparison similar to that of figure 4(a) with profiles con- 
verted to the x,y, z system using data from reference 12. Equations (29) 
do not afford a good fit in figure 4(b), although the equations do 
represent the general behavior of the velocities. This comparison is 
further illustrated in figure 5, where the value of g(y/b) as obtained 
by converting the profiles of reference 12 to the x,y,z system at indi- 
cated points is compared with (l-y/b) 2 . 


With the use of relations (29), H, K, J, and L are computed 
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H « 
J = 
K = 
L = 

which for n = 7 give 


2+n 

n 

n 2 (lln+7) 


{ 2n+l) ( 3n+l) [ 3n+2 ) 

\ 

2n^(2+n) , 
(2n+l)'(3n + l) 

6n 4 

/ 

( 3n+2) ( 2n+l) ( 5n+2) 

J 

H = 1.2857 
J a 0.5423 


K = 2.6727 


L = 1.1285 

J 


(30) 


(30a) 


Averaging the values along line III of Gruschwitz data (fig. 2) 
results in 


H = 1.37 


J a 0.550 
K = 2.43 


> (30h) 


L a 0.968 


J 


This relatively good agreement between the two sets of values is 
also indicative of the over -all fit of the assumed expressions for G 
and g to the data of reference 9. 


The additional relation that is needed for, the solution of the 
momentum-integral equation is the expression for surface friction. In 
reference 9, Gruschwitz demonstrates that von K^rmdn's friction law 

1 

, ^ . 0.01255 (5^)* (31) 

appears to be valid in the three-dimensional boundary layer as well. 
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Substituting relations (26), (28), and (31) in equations (23) . gives 

^ + Ti< 2+H > +j<d0v ^ 


+ Tfl 4 °V Tcfl _ 'OjX 

3z~ + J0 x 3^ ‘ 4 u" Je0 x “ — = ~ 


(32a) 


pU 


and 


2 L 0 X € gf + L€ ^ + (K_J)0 X ^ + ( K - J ) £ + gjjr + 

T 5x ' ’ ~ Tr Ux = -t — ( 32c ) 


ao. 


Be 


be. 


3 x au 


Oy. ^TT O 2l 

2(K-J)€ JE - (1+L€ 2 ) 


U 


pu 


Because of tbe form of the relation for friction, an advantageous 
transformation of variable is ^ 1 

0 = 0 x R 4 (0 x ) 

in order to "remove" the Reynolds number from the equations. 

With the use of equations (33) and (12), two nonlinear partial 
' differential equations for 0 and € applying along z = 0 are 
obtained from equations (32), 


(33) 


4 a© . 4 . T a© 

5S 5 J£ 


+ J©|| + 




au ■ j „ 

3x + 5 c€ 


^ J€ 0 -0.01255|= C 


(34a) 


and 


I (K - J)e i + 1 L * 2 i + (K - J) ® s + sf + 

[§ (K-J)e S - (| Lc 2 -1-Hjc0 - (f^-s) + 0.01255c] = 0 (34c) 

As shown in the appendix, these equations can be either hyperbolic, 
parabolic, or elliptic, depending on the shape of .the velocity profiles 
existing in the boundary layer. . For u = U(y/s). 1 / n and g = (l-y/&) , 
the equations are always elliptic • 
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Simplification of equations and approximate solution. - The rela- 
tive importance of the various derivatives in equations (34) is 'now- 
determined. First, 0 and € are assumed to he quantities that are 
.smaller than- one, which can be accomplished simply by referring all 
lengths to a total path length S and a to 45°. As a result of this 
assumption, all derivatives in 0 and € become of the order of mag- 
nitude of 0 or 6. Equation (34a) . is then divided by 4/5, so that the 
coefficient of b€)/bx is 1 and the coefficients of d©/dz and bc/bz 

are J€ and ^ J0 that is, of the order of magnitude of € or 0, 


respectively. In a similar way equation (34c) is divided by (K-J)0 
so that the coefficients of dQ/bx and bc/bx become — c/3 and 1 

(order of magnitude of one), respectively, and those of <3©/dz and 
bi/bz become — L and 2Lc (order of magnitude of e), respect- 
ively. Then, if e is small as compared with tan 45° and 0 is small 
when compared with S, all terms of the order of magnitude of c 2 , 0e, 
or 02 may be neglected, which gives 


b& 

3£ + 


5H+9\ 1 bH 
K 4 / U 3x 


0 = 0.01569 


(35a) 


and 


Sc 

3x + 


S0 9 1 Su 0.01255 1 
cE + 5 U 3x + (K-J) 0 


"T^.[ (1+B)c4H Tr] ^ 


These two expressions show that the primary changes in 0 and € in 
equations (34a) and (34c) occur only in the x-direction and thus the 
description of the phenomena only at z = 0 is justified. 

A solution of equations (35) can now be obtained by successive 
approximations because U, ojy, and c are assumed to be known func- 
tions of x. First, equation (35a) is solved. 



The values of 0-j-(x) are then used in an approximate solution of - equa- 
tion (35c) 
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4.9 
5 . ,5 


. \ vy 

(®l\ U i\ 1 1 

€ I -€ i\©i/ \ U/ 1^7 £9 

( J-K) (0 i ) 5 U 5 E 1 (x) 



r n 4 9 

f • ^ 5 5 

p+HjctHyj^U E^( 


x) dx 


where 


E x (x) = 


0.01255 

e -IWj xi ei 



(37) 


(38) 


With Cj(x) known, the approximation for © could in turn, he 

improved by again solving equation (35a). If a £rows large along the 
path, however, it is more advantageous to consider the following 
equation: 


a® + J€ a© + 

+ J€ 33 + 


l 5H+9 \ 1 
V 4 / U 


3U J 18 “y 

33 + I c€ " T Je 1 T 


© = 0.01569 


(39) 


It is thus hoped that negleetion of cte/dz will not affect the 
accuracy of the solution to a very marked extent. The solution of equa- 
tion (39), which may he obtained by the method of Lagrange, is 

''x /5H+9\ 

'x ' 


0 

II 


where 




E 3 (x) 


( 5H+9 \ + ( 5H+9 \ 

A’ 4 /Eo(x) lA 4 / 


0.01569 E 3 (x) « J IT 4 /E^x) 


- 1 uV" 4 /E 2 (x) U\ 4 /E 2 (x) J X i 


E 3 (x) 


E 2 (x) — e 



dx 

(40) 

(41) 


E,(x) * e 


T J l x . e! ? to 


(42) 
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and 



is an arbitrary function satisfying the boundary condition; when x = x^ 
and z = 0, then 0jj = 0.^. Setting 



at x = x^ gives A = 0. In addition, for z = 0, A is of the order 
of magnitude of £. 

Expansion of \jr in Maclaurin’s series about A(x^) yields 

A 2 

♦ (A) = t(0) + A*'(0) + 27 r '(0) + . . . 


Because there is only one boundary condition, it is possible to determine 
o nly one of the constants in this expansion; consequently, \jr cannot be 
uniquely established. The fact that A is of the order of magnitude of 
e, however, suggests that the assumptions made for ^'*(0), i|/'''(0), 

and so forth, are successively less important. Thus these derivatives 
may arbitrarily be expressed by a single constant, 

( 5H+9 ^ / 2.2 \ ( 5H+9 \ 

t(A) - %U i V 4 J {l+A\ + ^27- + • • -) - ®i u i 4 ' / e AA ( 4 3 ) 

where A from purely dimensional considerations must have the dimen- 
sions of 2"^-. From expressions (41) and (42), it is suspected that 

; < 44 > 

where B must be obtained from the experimental measurements. 

It should be noted that because i|r(A) cannot be uniquely determined 
other functions of A satisfying the single boundary condition could be 
used as well. The function e^X i s chosen only because it is convenient 
to use and parallels the expressions (41) and (42). This arbitrariness 
of the functional form of \|r and the value of A is due to the con- 
sequences of ass uming z = 0, and thus it is probably not advisable to 
carry any further approximation for C-j-j, and so forth. 
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In solving equations (36) to (44), either set of values for H, K, 
J , and L may be used. Because the averaged values (30b) are probably 
more representative, having been obtained by evaluating experimental 
data at a number of different positions, it is advantageous to use 
these values in computations. 


COMPARISON WITH EXPERIMENT 

In order to check the validity of the approximate solution, the 
boundary layer along four streamlines of reference 9 was computed and 
compared with the measured values. The designation of the streamlines 
and data points is illustrated in figure 2. Because the data were 
taken along curves I to V of figure 2, the computation along a stream- 
line requires first an interpolation among the various data points. As 
a result of this interpolation, the computations could not be carried 
through the full length of each streamline. Values (30b) were used for 
quantities H, K, J, and L. The constant B was obtained by fitting 
along streamline "B" the solution for 0 -q, so that at x = S, 

©II'-© measured. In this manner, the value of B was found to be 

38.5. This value was then used in computations of streamlines "A", 

"C", and "D". It is noted that B = 38.5 ~ 7R 1 / 4 (e x ) i , although 
justification for such a dependence cannot be made. In all integrations 
Simpson l s rule was used. 

The results of the computations are plotted in a nondimens ional 
form and compared with Interpolated measured values in figures 6 and 7. 

A study of these figures reveals a fair quantitative agreement between 
the measured and estimated values of 0 and a. As the values of 
€— ^tan 45° (fig. 7) the first approximation for 0 in figure 6 becomes 
progressively worse, which is remedied by the second approximation. The 
poorest agreement is obtained along streamlines "A" and "D", which 
because they are closest to the walls , might be affected by the flow in 
the corners of the duct. Streamline "D" especially may thus be 
affected inasmuch as Gruschwitz mentions the existence of separation 
on the convex wall. 

The fair quantitative agreement with the measured values is not to 
be interpreted as a conclusive check of. the validity of the procedure 
and the assumed values in all cases of three-dimensional boundary -layer 
flow. The suggested procedure simply represents the best that can be 
done in view of the meagerness of the available data. Because the 
Gruschwitz data do not involve the effects of uniform angular velocity 

and because the variations in ^ and jj are small, it could 

be maintained that this check of the procedure has been carried out 
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on a somewhat special case. For that reason, it is desirable that 
additional experiments be carried out in setups that eliminate, the 
present shortcomings. Larger variation of Reynolds number should 
also be used. With additional experiments, a modification of the 
values of H, K, J, L, and B, together perhaps with some 
refinements of the procedure will be in order. It might be well 
to remember at such time, that because of the necessary empiricism 
involved (which results from the very limited knowledge of turbulent 
phenomena), long and tedious computations would rarely be worthwhile. 


CONCLUSIONS 

The following conclusions can be drawn from an analysis of the 
three-dimensional momentum- integral equations and a comparison of the 
numerical results with the Gruschwitz data for turbulent boundary layer: 

1. Within the boundary layer the static pressure can vary at most 
by an amount of the order of magnitude of the boundary-layer 
thickness 5 . 

2. It is possible to generalize the velocities in the boundary 
layer by use of two characteristic quantities 8 and € where e is 
the tangent of the angle enclosed by the direction of the resultant 
skin-friction stress and the direction of the flow outside the boundary 
layer . 

3. When the generalized boundary-layer momentum-loss thickness 0 
is small as compared with the total path length and € is small as 
compared with tan 45°, the primary changes in 0 and € occur along 
the streamline of the flow outside the boundary layer. 

4. The three-dimensional boundary-layer momentum- integral equations 
can be either hyperbolic, parabolic, or elliptic, depending on the rel- 
ative magnitude of the parameter MN, which in turn depends on the 
shape of the velocity profiles existing in the boundary . layer . The 
power-law profile when used with the correction function g = (1-y/fe) 2 
always results in elliptic’ equations. 

5. The approximate solution of the three-dimensional momentum- 
integral equations shows a fair quantitative agreement with the values 
measured by Gruschwitz. 

6. Additional experimental data are necessary to establish more 
generally applicable values for form parameters H, K, J, and L 
and B, the constant used in the second approximation for 0. 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 

Cleveland, Ohio, November 1, 1950. 
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APPENDIX - ADDITIONAL REMARKS ON THREE-DIMENSIONAL BOUNDARY -LAYER 
MOMENTUM-INTEGRAL EQUATIONS' 


In order to obtain the approximate* solution of equation (34), it 
was shown by comparing the relative order of magnitude of the coeffi- 
cients that some of the terms may be neglected. Care must be taken with 
such simplifications inasmuch as various implications of the equations 
in question may be obscured by this procedure. For this reason, aside 
from the approximate solution, the character of equations (34) was also 
investigated in detail. 

By use of the procedure outlined in reference 18 (p. 38) , along 
z = 0 the system of equations (34) is found to be hyperbolic when 
MN >1, elliptic when MN < 1, and parabolic when MN = 1. 


where 


* s Tic^yj 


0 XZ 0 zx 


(45) 


Because L, K, and J are functions of G and g, the char- 
acter of equation (34) depends on the .shape of the velocity profiles in 
the boundary layer. 

It should be noted that when MN = 0, then L = 0, which is only 
possible if G « 0 or g = 0, and in turn u = 0 or w = 0. If the 
trivial case u = 0 is neglected, it is established that when w = 0, 

C = 0 as well. But for € = 0 and w = 0, equations (34) reduce to a 
special case ' 

I g + ! tW) ■ °- 01255 (46a) 

and 

2»jr 1+H _ _ 0.01255 , (46c) 

U H ® . 

Here equation (46a) is an ordinary two-dimensional boundary-layer 
momentum-integral equation for 0 and equation (46c) is a relation that 
evidently must exist among U, du/dz, ay, and 0, when € = 0 and 
w = 0. 
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When the equations are elliptic, no real characteristic direction 
can he found. When only one characteristic direction exists, the equa- 
tions are parabolic and in the hyperbolic case two characteristic 
directions through each point of the xz -plane are obtained. For the 
parabolic case then, 



Je 


• and for the hyperbolic case, 


dz _ f L J_*j L 2 -J(K-J)L 

dx =+ “ (K-J) 


dz , „ L - L 2 -J.(K-J)L 

to = £-■“ c (k^jJ 


and the characteristic lines are asymetric with respect to the x-axis. 

In order to determine whether elliptic, parabolic, or hyperbolic equa- 
tions apply, the magnitude of MET is computed. Substituting from 
expressions (30), MN ' is obtained in terms of n: 

6(3n+l)(3n+2) 1 108n+96 

" (Sn+2) (lln+7) " ' 55 ' 55 ( 5 5 n 2 + 57 n+14 ) 

This equation shows MN to be a monotonically increasing function of n. 
For n = 0, ' 

6 

MN - y 

and 


lim MN = 55 

n -4 oo 

These results indicate that a so-called power -law profile when used 
with g - (l-y/ 8) 2 always results in equations that although 
elliptic are very near to being parabolic. Using values (30b), 

MN = 0.936 

which again indicates an elliptic character of the equations. It should 
be remembered, however, that the assumption for G and g were made 
on the basis of only one set of data; consequently there is no assurance 
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that the velocity distributions existing in the boundary layer will 
always give the same values of MR. In fact, it is generally more 
likely that they will not give the same values of MN. Some indication 
of the variation of MR may already be obtained from f igure 8, where 
its value was plotted for each Gruschwitz data point. At points 5, 6, 
and 7, however, small values of w may have prevented an accurate 
determination of L and as a consequence MN = 0 there. The value of 
MN in figure 8, varies within the limits 0.65<MN< 1.2, with the bulk 
of the points indicating that MN-0.95. 

On the basis of the preceding discussion, there is some evidence 
of the equations being parabolic, elliptic, and hyperbolic in the tur- 
bulent boundary layer. It is interesting to note that generally (as 
in supersonic and subsonic flow, for instance), these hyperbolic and 
elliptic regions have their counterpart in physical phenomena. Thus 
some essential differences might exist in the process of momentum 
transfer between the hyperbolic and elliptic regions. These differ-, 
ences cannot now be ascertained because' first equations similar to (34) 
with z / 0 would have to be obtained, and there is no mention of any 
irregularities in the behavior of the flow in reference 9. When 
additional experiments are made, however, it would seem advisable to 
closely study these two mathematical regions in order to obtain some 
indication of the physical make-up of their differences. 
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X = / cos p dx + z sin P 
J 0 

Y=y,Z=Z^+z cos P 


Z^ = constant 
dp 

c = dx- 


■r 


sin P dx 


“0 

P = P(x) 


Figure 1. - Transformation from Cartesian coordinates 
X.Y.Z to orthogonal curvilinear coordinates x,y,z. 




NACA TW 2310 



(a) Channel and measuring plate seen from "below. 


18 19 20 



(h) Measuring plate seen from "below, showing 
point and streamline designations. Shaded 
sections indicate regions of pot ential- flow 
"breakdown. 

Figure 2. - Experimental set-up of Gruschwitz 
(from figs: 1 and 5 of reference 9). 









Assumed profiles 
Data of reference 12, position 
converted to x,y,z system 
Data of reference 12, position 
converted to x,y,z system 
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Figure 4. - Concluded. Comparison of assumed G and Gg with experimental Telocity profiles. 

Assumed profiles: G = (y/6) 1 ^ 7 ; Gg = (y/6)^ (l _ y/S) • 




Generalized momentum thickness , ©/0 



Generalized momentum thickness, 0/0^ 



(b) Streamline "B"j S = 34.75 inches; 0 ± « 0.289 inch. 

Figure 6. - Continued. Comparison of calculated and measured 
generalized momentum thickness. Experimental data from 
reference 9. 




Generalized momentum thickness , 0/0. 
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Calculated, approximation I - 
Calculated, approximation II 


Distance along streamline , x/S 
(d) streamline "D"; S = 21.25 inches; 0^ =0.244 inch. 

Figure J5. - Concluded. Comparison of calculated and measured 
generalized momentum thickness . Experimental data from 
reference 9. 
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Figure 7. - Concluded. Comparison of calculated and measured 
"boundary -layer deflection at vail. Experimental data from 
reference 9. 
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Position of data points 


Figure 8. - Values of parameter MN for data of reference 9. 


NACA-Langley - 3-12-51 - 1025 



